Introduction to the Bethe Ansatz II 
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I. INTRODUCTION 

Quantum spin chains are physically realized in quasi- 
one-dimensional (ID) magnetic insulators. In such ma- 
terials, the magnetic ions (for example, Cu ++ and Co ++ 
with effective spin s = \, Ni ++ with s = 1, and Mn ++ 
with s = | ) are arranged in coupled linear arrays which 
are well separated and magnetically screened from each 
other by large non-magnetic molecules. Over the years, 
magneto-chemists have refined the art of designing and 
growing crystals of quasi- ID magnetic materials to the 
point where physical realizations of many a theorist's pet 
model can now be custom-made. 

Most prominent among such systems are realizations 
of the ID spin- 1/2 Heisenberg antiferromagnet, 



N 



Ha — J ^ S n • S n 



(1) 



n=l 



the model system which inspired Bethe to formulate the 
now celebrated method for calculating eigenvalues, eigen- 
vectors, and a host of physical properties. Interest in 
quasi- ID magnetic materials has stimulated theoretical 
work on quantum spin chains from the early sixties un- 
til the present. Some of the advances achieved via the 
Bethe ansatz emerged in direct response to experimen- 
tal data which had remained unexplained by standard 
approximations used injpany-body theory. 

In Part I of this seriesEl we introduced the Bethe ansatz 
in the spirit of Bethe's original 1931 paper.cl The fo- 
cus there was on the Hamiltonian Hp = —Ha, where 
the negative sign makes the exchange coupling ferromag- 
netic. The ferromagnetic ground state has all the spins 
aligned and is (JV+l)-fold degenerate. The reduced rota- 
tional symmetry of every ground state vector relative to 
that of the Hamiltonian reflects the presence of ferromag- 
netic long-range order at zero temperature. One vector 
of the ferromagnetic ground state, \F) = || ■ • • T), in the 
notation of Part I, serves as the reference state (vacuum) 
of the Bethe ansatz. All other eigenstates are generated 
from \F) through multiple magnon excitations. 

We also investigated some low-lying excitations which 
involve only a small number of magnons. For example, 
the Bethe ansatz enabled us to study the properties of a 
branch of two-magnon bound-states and a continuum of 
two-magnon scattering states, and gave us a perfect tool 
for understanding the composite nature of these states in 
relation to the elementary particles (magnons). 



In this column we turn our attention to the Heisen- 
berg antiferromagnet Ha- All the eigenvectors remain 
the same as in Hp, but the energy eigenvalues have the 
opposite sign. Therefore, the physical properties are very 
different, and the state \F) now has the highest energy. 
Our immediate goals are to find the exact ground state 
\A) of Ha, to investigate how the state \A) gradually 
transforms to the state \F) in the presence of a magnetic 
field of increasing strength, and to prepare the ground 
work for a systematic study of the excitation spectrum 
relative to the state \A). 

As in Part I we will emphasize computational aspects 
of the Bethe ansatz. The numerical study of finite sys- 
tems via the Bethe ansatz is akin to a simulation in many 
respects, and yields important insights into the underly- 
ing physics. 



II. GROUND STATE 

What is the nature of the ground state state \A)1 How 
do we find it? Is its structure as simple as that of \F)1 
An obvious candidate for \A) is the Neel state. The two 
vectors 



|M) = |UT •••!>, |AA 2 ) = |m---T) 



(2) 



reflect antiferromagnetic long-range order in its purest 
form just as the vector \F) does for ferromagnetic long- 
range order. Henceforth we assume that the number of 
spins N in (|l|) is even and that periodic boundary condi- 
tions are imposed. 

Inspection shows that neither IA/2), nor the trans- 
lationally invariant linear combinations, \Af±) = (|A/i) ± 
IA/2})/ v^, are eigenvectors of Ha- In the energy expec- 
tation value (Ha), the Neel state minimizes (S^S^+i) 
but not {S£S% +1 ) and {S»S% +1 ) (Problem 0). A state 
with the full rotational symmetry of Ha can have a 
lower energy. Like |A/±), the ground state \A) will be 
found in the subspace with Sf, = ^2 n S* = 0. Because 
of their simplicity, the Neel states are convenient start- 
ing vectors for the computation of the finite- N ground 
state energy and wave function via standard itera- 
tive procedures (steepest-descent and conjugate-gradient 
methods) B However, here we take a different route. 

In the framework of the Bethe ansatz, all eigenstates of 
Ha with Sj. = can be obtained from the reference state 
\F) by exciting r = N/2 magnons with momenta fe, and 
(negative) energies — J(l— cos hi). The exact prescription 
was stated in Part I. Each eigenstate is specified by a 



1 



different set of N/2 Bethe quantum numbers {A^}. The 
momenta ki and the phase angles Qij in the coefficients 
(128) of the Bethe wave function (I2Z) result from the 
Bethe ansatz equations (133) and (135). □ A finite- N study 
indicates that the ground state \A) has real momenta ki 
and Bethe quantum numbers (Problem |^) 

{A i }A = {l,3,5,...,JV-l}. (3) 

In Part I we worked directly with ki and fty, which 
represent physical properties of the elementary particles 
(magnons) created from the vacuum \F). At the opposite 
end of the spectrum, the state \A), generated from \F) 
via multiple magnon excitations, will be configured as a 
new physical vacuum for Ha- The entire spectrum of 
Ha can then be explored through multiple excitations of 
a different kind of elementary particle called the spinon. 

Computational convenience suggests that we replace 
the two sets of variables {ki} and {Oij} by a single set 
of (generally complex) variables {zi}. If we relate every 
magnon momentum ki to a new variable Zj by 

ki = TT-(p(Zi), (4) 

via the function 4>{z) = 2 arctanz, then the relation (133) 
between every phase angle Oij for a magnon pair and the 
difference Zi — Zj involves 4>{z) again: 

i:j = 7rsgn[3?(zi - Zj)] - 4>[{zi - Zj)/2]. (5) 

Here 3?(x) denotes the real part of x, and sgn(y) = ±1 
denotes the sign of y. Relations (|J) and (|J) substituted 
into (I35-) yield the Bethe ansatz equations for the vari- 
ables zfE 

N<l>{zi)=2nIi + ^if>[{zi-z j )/2], i = l,...,r. (6) 

The new Bethe quantum numbers F assume integer val- 
ues for odd r and half-integer values for even r. The 
relation between the sets {Xi} and {/j} is subtle. It de- 
pends, via sgn[3?(zi — Zj)], on the configuration of the 
solution {z^ in the complex plane. For the ground state 
| A), we obtain (Problem ||a) 




Given the solution {z%, . . . , z r } of Eqs. (|J) for a state 
specified by . . . , I r }, its energy and wave number are 
(Problem |b) 

r 

{E-E F )/J = ^e{zi), (8a) 

i=l 

*=i>-0(*)] =^-^$>, ( 8b ) 

i=l i=l 

with Ep — JN/A. The quantity 4>(zi) is called the 
magnon bare momentum, and e{z/) = dki/dzi — —2/(1 + 



zf) is the magnon bare energy. The Bethe wave function 
(127) is obtained from the {z^ via (@) and (|). 

The ground state \A) belongs to a class of eigenstates 
which are characterized by real solutions {zi} of the 
Bethe ansatz equations. To find them numerically we 
convert Eqs. (Q) into an iterative process: 

4 n+1) = **(~i< + ± E^[^ (n) - «S B) )/2]) • ( 9 ) 

Starting from z^ -* = 0, the first iteration yields z^ — 
tan(7r Ii/N). Convergence toward the roots of (||) is fast. 
High-precision solutions {zi} can be obtained on a per- 
sonal computer within seconds for systems with up to 
-/V = 256 sites and within minutes for much larger sys- 
tems (Proble m |4| ). The ground-state energy per site in- 
ferred from (paf) is listed in Table I for several lattice 
sizes. 

TABLE I. Numerical results for the energy per site of the 
ground state \A) relative to the energy Ef/JN = 1 of the 
state \F) for various values of iV obtained via n max iterations 
of (^|). The CPU time quoted is for a Pentium 130 computer 
running GNU C on Linux. The exact result for N — * oo is 
-In 2. 



N 


(E A - E F )/JN 


^max 


CPU-time [sec] 


16 


-0.696393522538549 


38 


0.01 


64 


-0.693348459146139 


83 


0.08 


256 


-0.693159743366446 


195 


2.38 


1024 


-0.693147965376242 


483 


104 


4096 


-0.693147229600349 


1256 


4695 


oo 


-0.693147180559945 







In preparation of the analytical calculation which pro- 
duces (Ea — Ep)/JN for N — > oo, we inspect the Zi- 
configuration for the finite- N ground state \A) obtained 
numerically. All roots are real, and their values are sorted 
in order of the associated Bethe quantum numbers 7j. 
The 16 circles in Fig. |l|(a) show z, plotted versus Ii/N 
for N = 32. The solid line connects the corresponding 
data for N = 2048. The line-up of the finite- TV data 
along a smooth monotonic curve suggests that the solu- 
tions of (|^) can be described by a continuous distribution 
of roots for N — > oo. 

We give the inverse of the discrete function shown in 
Fig. |l](a) the name Zj^(zi) = F/N and rewrite Eqs. (|^) 
in the form: 

2nZ N ( Zi ) = 4>{ Zi ) - i <f>[{zi - Zj)/2] . (10) 

For N — > oo, Zff{zi) becomes a continuous function Z{z) 
whose derivative, gq{z) = dZ/dz, represents the distri- 
bution of roots. In Eq. (JToJ) the sum (1/N) J^ . . . . is 
replaced by the integral f dz'crp (z') . . . Upon differenti- 
ation the continuous version of (nm becomes the linear 
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(a) 



- 



-0.2 





Z N =I/N 



0.2 




Z N =I,/N 



FIG. 1. Ground state \A). (a) Solutions of Eq. (|) and (b) 
magnon momenta in Eq. (Q) plotted versus the rescaled Bethe 
quantum numbers (R) for A' = 32 (circles) and TV = 2048 
(lines) . 



integral equation 

2tt<j (z) = -e(z) - (K * a )(z) 



(11) 



with the kernel K(z) = 4/(4 + z 2 ); (K * a Q )(z) is short- 
hand for the convolution f_°° dz' K(z— z')ao(z'). Fourier 
transforming Eq. ( pd| ) yields an algebraic equation for 
(To(u) = dz e luz <To(z). Applying the inverse Fourier 
transform to its solution yields the result 



<t (z) = -sech(7rz/2). 



(12) 



The (asymptotic) ground-state energy per site as inferred 
from (Ka 



Ea-E f 
JN 



-t-oo 



dz e(z) <jq(z) 



In 2, 



(13) 



is significantly lower than the energy expectation value 
of the Neel states (Problem [j]). The state \A) has to- 
tal spin St — (singlet). Unlike \F), it retains the full 
rotational symmetry of. (Q) and does not exhibit mag- 
netic long-range order .0 The wave number of \A), ob- 
tained from (|l]), is Ua = for even N/2 and fc^ = 7r 
for odd-A^/2. The important result ( |T^ ) was derived by 
HulthcnQ in the early years of the Bethe ansatz. 

In Fig. 0(b) we plot the magnon momenta fcj of \A) as 
inferred from ([|) versus U/N for the same data as used 
in Fig. |l|(a). The smoothness of the curve reflects the 
fact that the state \A) for N — > oo can be described by a 
continuous /^-distribution (Problem |J) 



Po(k) = 









8 sin 2 — cosh 




cot-j 


2 * 







1 -1 



(14) 



III. MAGNETIC FIELD 

In the presence of a magnetic field h, the Hamiltonian 
(|l|) must be supplemented by a Zeeman energy: 



H = Ha — hS^ 



(15) 



The two parts of H are in competition. Spin alignment 
in the positive z-direction is energetically favored by the 
Zeeman term, but any aligned nearest-neighbor pair costs 
exchange energy. Given that [Ha, S t ] = 0, the eigenvec- 
tors are independent of h. The 2St + 1 components (with 
\S T \ < St) of any Sr-multiplet fan out symmetrically 
about the original level position and depend linearly on 
h. 

The largest downward energy shift in each multiplet is 
experienced by the state with S T = St, and that shift 
is proportional to St- The state \A), which has St = 0, 
does not move at all, whereas the state \F) with St = 
N/2 descends more rapidly than any other state. Even 
though \F) starts out at the top of the spectrum, it is 
certain to become the ground state in a sufficiently strong 
field. The saturation field h$ marks the value of h where 
\F) overtakes its closest competitor in the race of levels 
down the energy axis. 

The pattern in which levels with increasing S T become 
the ground state of H as h increases depends on their 
relative starting position along the energy axis. From 
the zero-field energies of this set of states, we will now 
determine the magnetization m z = S T /N of the ground 
state as a function of ha 

The Bethe quantum numbers of the Lowest state with 
quantum number S T = N/2 — r > ared 



l + 2i- 



N 
~2 



(16) 



as can be confirmed by finite- TV studies of all states in the 
invariant subspaces with r = 1, . . . , N/2. Using the iter- 
ative process (|9|) , we can determine the energies of these 
states with high precision (Problem ^J). The red circles 
in Fig. ||(a) represent the quantity [E(S T ) — Ef]/JN for 
N = 32. The solid line connects the corresponding re- 
sults for N = 2048. For the finite- TV analysis it is impor- 
tant to note that both the level energies E{S T ) and the 
level spacings E(S T ) — E(S T — 1) increase monotonically 
with S T . 

At h ^ 0, all of these levels experience a downward 
shift of magnitude hS T . All spacings between adjacent 
levels shrink by the same amount h. The first level 
crossing occurs between the state \A) with S T = and 
the state with S T = 1, which thereby becomes the new 
ground state. Next, this state is overtaken and replaced 
as the ground state by the state with S T = 2 and so forth. 
The last of exactly N/2 replacements involves the state 
with S T = N/2 - 1 and the state \F) with S T = N/2. 
Their energy difference in zero field is 2J independent 
of N (see Part I). Consequently, the saturation field is 
h s = 2 J. 
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FIG. 2. (a) Energy of the lowest states with a given Sf> at 
h = 0. (b) Magnetization in the ground state of H. 

The magnetization m z grows in N/2 steps of width 
l/N between h = and h = hs- In Fig. ^(b) we plot 
m z versus h for two system sizes based on data obtained 
numerically. The blue staircase represents the results for 
N = 32. For N = 2048 the step size has shrunk to well 
within the thickness of the smooth curve shown. 

An important observation is that the midpoints of the 
vertical and horizontal steps of the m z (h) staircase (red 
dots) lie very close to the limiting curve. This behavior 
made it possible to extract quite accurate magnetization 
curves for various spin chain models from very limited 
datafl 

Different scenarios are conceivable. For example, if the 
levels were arranged in the same sequence as in Fig. |^(a) 
but with spacings increasing from top to bottom, then 
the state with S T = would be replaced directly by the 
state with S T = N/2. 

In the thermodynamic limit, the energy per site of the 
lowest level with given S T becomes the internal energy 
density at zero temperature, 



u{m z ) 



lim 



E\Sj<) — Ep 

In ' 



(17) 



From ( |17| ) we obtain, via the thermodynamic relations, 
du dm z ( d 2 u x ' 



dm. 



Xz 



dh 



\dm 2 z 



(18) 



The function m z (h) shown in Fig. ||(b) is the inverse of 
the derivative of the function u{m z ) plotted in Fig. |](a). 
The slope of m z (h) represents the longitudinal suscepti- 
bility Xzz- In a finite system, where m z = S T /N varies 
in steps of size l/N, Eqs. Jig ) are replaced by 



h(m z ) = E(S T ) - E{S T - 1), 



Xzz(m z ) 



l/N 



E(S y T + 1) - 2E{S T ) + E(S T - 1) ' 



(19a) 
(19b) 



The data in Fig. |^(b) indicate that Xzz(h) has a 
nonzero value at h = 0, grows monotonically with h, 



and finally 
initial value 



fir 



rges at the saturation field h = hs- The 



Xz 



7T J 



(20) 



turns out to be elusive to a straightforward slope analysis 
because of a logarithmic singularity which produces an 
infinite curvature in m z (h) at h = (Problem 0). The 
divergence of Xzz(h) at hs is of the type (Problem ||) 



Xzz(h) 



h^hs 1 



1 



27T 



h) 



(21) 



The characteristic upwardly bent magnetization curve 
with infinite slope at the saturation field is a quantum 
effect unreproducible by any simple and meaningful clas- 
sical model system. The Hamiltonian (Q), reinterpreted 
as the energy function for coupled three-component vec- 
tors, predicts a function m z (h) which increases linearly 
from zero all the way to the saturation field. 



IV. TWO-SPINON EXCITATIONS 

Returning to zero magnetic field, let us explore the 
spectrum of the low-lying excitations. From here on, the 
ground state \A) (with = St = 0) replaces \F) (with 
S T = St = N/2) as the new reference state for all ex- 
cited states. The Bethe quantum numbers (0), which 
characterize \A), describe a perfectly regular array on the 
/-axis as illustrated in the first row of Fig. |[ This array 
will be interpreted as a physical vacuum. The spectrum 
of Ha can then be generated systematically in terms of 
the fundamental excitations characterized by elementary 
modifications of this vacuum array. 
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FIG. 3. Configurations of Bethe quantum numbers h for 
the JV = 32 ground state (top row) and for one representative 
of three sets of two-spinon excitations with energy E < ^ t ' St \ 
Each gap in the /^-configurations of rows two and three (green 
full circles) represents a spinon. Each gap in row four (green 
open circles) represents half a spinon. The blue circle rep- 
resents the Bethe quantum number associated with a pair 
of complex conjugate solutions, whereas all black circles are 
associated with real solutions. 

In the subspace with S T — 1, a two-parameter set 
of states is obtained by removing one magnon from the 



4 



state \ A). In doing so we eliminate one of the N/2 Bethe 
quantum numbers from the set in the first row of Fig. ^ 
and rearrange the remaining Ii in all configurations over 
the expanded range —\N < Ii < iJV. Changing by 
one means that the Ii switch from half-integers to inte- 
gers or vice versa. The number of distinct configurations 
with Jj+i — Ii > 1 is ^N(N + 2). A generic configuration 
consists of three clusters with two gaps between them as 
shown in the second row of Fig. ||. 

The two gaps in the otherwise uniform distribution 
of Bethe quantum numbers can be interpreted as ele- 
mentary particles (spinons) excited from the new phys- 
ical vacuum. The position of the gaps between the Jj- 
clusters determine the momenta ki , of the two spinons, 
which, in turn, add up to the wave number of the two- 
spinon state relative to the wave number of the vacuum: 
q = k - h,A = ki + ki- 

What remains to be done for the finite- N analysis is 
straightforward. We solve the Bethe ansatz equations 
via (H) for all the ^-configurations just specified. A plot 
of the two-spinon energies E — Ea versus wave number 
k — kA for N = 16 as inferred from the solutions {z^ via 
(paf) and (8b) is shown in Fig. 



(red circles) . Also shown 
are the Bethe quantum numbers for each excitation. The 
pattern is readily recognized and extended (via reflection 
Ii <-> —Ii) to the other half of the Brillouin zone. We have 
drawn (in blue) the corresponding two-spinon excitations 
for N = 256. The jN(±N + 1) = 4160 dots in the range 
< q < 7r produce a density plot for the two-spinon 
continuum which emerges in the limit N — > oo. 



of the two-spinon continuum,t2l 

7T 

e L{q) = - J\sinq\, eu{q) 



7rJ|sin||, (22) 



starts from the Bethe ansatz equations (10). We set 
r = N/2 — 1 and use the Bethe quantum numbers from 
the second row of Fig. |[ When we replace the sum by 
an integral, we must account for the gaps between the 
ii-clusters by two 1/iV-corrections. The result (after dif- 
ferentiation) is an integral equation which differs from 
( |ll| ) by two extra terms related to the I^-gaps: 

s(z)-(K*a)(z) + ^^2K(z-z l ). (23) 



2tt(t(z) 



1=1,2 



If we write a(z) = (Jq{z) + &i{z) + 02(2:), where <Jq{z) is 
the solution ( |i~2| ) of Eq. ( |Tl| ) , the two-spinon corrections 
are solutions of 

2na l (z) = ^K(z-zi)-(K *<n)(z), 1=1,2. (24) 



Because Eqs. (f24|) have the same integral kernel K as 
Eq. ([ll]), the solutions of all three equations can be ex- 
pressed by the same resolvent R, 



2irai(z)= gi {z)-{R* gi )(z), 1 = 0,1,2 



(25) 



where on(z) and gi(z), Q2(z) are the inhomogeneities of 
Eqs. (|ll| ) and (pif), respectively. The resolvent, which 
does not depend on the inhomogeneity, satisfies 




FIG. 4. Two-spinon (triplet) excitations with St = St = 1 
for iV = 16 (red circles) and N = 256 (blue dots). 
The configurations of (integer) Bethe quantum numbers 
(-4 < h < ... < h < 4) for iV = 16 are symbolized by 
the black circles. 

Compare this set of two-spinon scattering states with 
the set of two-magnon states plotted in Fig. 2 of Part I. 
There we found that pairs of magnons form a continuum 
of two-magnon scattering states and a branch of two- 
magnon bound states. Two-spinon bound states exist as 
well and will be discussed in a later column. 

A derivation of the exact lower and upper boundaries 



2nR{z) = K{z) -{K* R)(z). 



(26) 



Because gi{z) = (l/N)K(z — zi) for I = 1,2 makes (|24|) 
equivalent to (|2^), we can write <Ji(z) = (l/N)R(z 
This is all we need to know about the solutions. 

Using Eq. ( pa] ) for the two-spinon energies, we must 
again correct for the two J^-gaps when we convert the 
sum into an integral: 



E — Ef — NJ 



dz e(z)a(z) — J e(zi). 

1=1,2 



(27) 



Subtracting Ea — Ep yields 
E-E A = -jJ2 ~ N (°i * £ )(°)] 

1=1,2 

= J i9o(zi) - (R*g )(zi)} 



1=1,2 



1=1,2 



~2 



5> ch (?) 

1=1,2 



(28) 



Now we must relate the J^-gaps, that is, the values z\, zi 
to the spinon momenta k\,ki- Starting from the config- 
uration in the first row of Fig. |^, we remove one Bethe 
quantum number at or near the center. The accompa- 
nying integer «-> half-integer switch of the remaining 
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opens a double gap at or near the center. Shifting sev- 
eral Ii from the left and right toward the center then 
produces the generic three clusters. This change in con- 
figuration applied to (pb|) yields the following two-spinon 
wave number relative to the vacuum: 



q = k-kA = ki + k 2 = TT-^ ^ Al, 



N/2-1 



N 



i=l 



- - 2n J dza (z) 



E 

1=1,2 

E [f ~ arctan ( sinh 



(29) 



or sinfc; = sech(7rz;/2), 1 = 1,2. These relations substi- 
tuted into Eq. (Eq) produce the two-spinon spectrum 



(30) 



E- E A = — J[j sinfcij + | sinfc 2 |] 
= 7rJ|sm(g/2)cos(p/2)| 



for < p < q and < q < 2~k — p. The spectrum is a 
two-parameter continuum with boundaries (p2|). 

Like magnons, spinoa*, carry a spin in addition to 
energy and momentum.liil Unlike magnons, which are 
spin-1 particles (see Part I), spinons are spin-^ parti- 
cles. In a chain with even N, where all eigenstates have 
integer- valued S T , spinons occur only in pairs. The spins 
si = ±4, / = 1, 2 of the two spinons in a two-spinon eigen- 
state of Ha can be combined in four different ways to 
form a triplet state (St = 1 5 S T — 0, ±1) or a singlet 
state (St — S^ — 0). They are described by distinct 
configurations of Bcthe quantum numbers. 

We have already analyzed the spectrum of the two- 
spinon triplet states with si = S2 = +| (St = 1, S T = 
+1), as specified by /^-configurations of the kind shown 
in the second row of Fig. ||. The two-spinon states with 
s\ = S2 = —1 (St = 1, S^ — —1) are obtained from 
these states by a simple spin-flip transformation. They 
have exactly the same energy-momentum relations. 

The remaining two sets of two-spinon states have si = 
— S2, that is, S T = 0. One of these sets contains triplets 
(St = 1) and the other singlets (St = 0). The former 
is obtained if we shift all Bethe quantum numbers in the 
second row one position to the right (7j — > 7j + |, i = 
1, . . . , r — 1) as shown in the third row and add one non- 
movable Bethe quantum number I r = |iV + | to the 
set (Problem [9]) . The number of distinct configurations 
is then the same as in row two and the corresponding 
states have the same wave number and energy. Symmetry 
requires that all three components (5f. = 0,±1) of the 
triplet states are degenerate. 

Whereas all two-spinon triplets are described by real 
solutions {zi}, the two-spinon singlet states (S T = St = 
0) are characterized by one pair of complex conjugate so- 
lutions z\ = u + iv, u — iv in addition to the real so- 
lutions Z3, . . . , zjv/ 2 E3 Finding the ii-configurations for a 
particular set of eigenstates with complex solutions (such 



as indicated in row four of Fig. ||) and then solving the as- 
sociated Bethe ansatz equations turns out to be a much 
more delicate task than was the identification and the 
calculation of purely real solutions. 

A straightforward adaptation of the iterative process 
(^) to complex Z{ will, in general, not converge toward a 
physically meaningful solution. Making the Bethe ansatz 
equations amenable to root-finding algorithms which are 
adequate for this task requires that we convert Eqs. (||) 
into a set of equations with real solutions. For the two- 
spinon singlets we arrive at the following set of N/2 equa- 
tions for {u, v, Z3, . . . , 2^/2} (Problem |l0| ): 
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, (31c) 



where ip(x) = 2atanh(x). Here we have replaced the 
first two Bethe quantum numbers I\ , I2 of Eq. (^) by a 
single (integer) Bethe quantum number l' 2 ' associated 
with z\ — z\. The Bethe cmantum numbers associated 
with real solutions, now renamed l\ \ , are constrained 

to the range (-jN - \ < < \N + \) as indicated 
in the fourth row of Fig. 0. This notation is used in 
preparation of a general classification of Bethe ansatz 
solutions (string hypothesis) which is valid for N — > 00 
and will be discussed in a later column. 

The dependence of the wave number on the Bethe 
quantum numbers for all states with one pair of complex 
solutions is then 
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(32) 



where r = N/2 for the two-spinon singlets considered 
here. 

Figure |^ depicts the energy (|8a| ) versus the wave num- 
ber ([52|) of the two-spinon singlets (red circles) for N = 
16 and < q < tt. We use the same representation as in 
Fig. [| and show the two-spinon triplets again for compar- 
ison. The number of two-spinon singlets identified here is 
smaller, namely i./V(iV— 2) over the entire Brillouin zone, 
than the number %N(N + 2) of two-spinon triplets iden- 
tified previously. Although the Ij-configurations of the 
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TABLE II. Two-spinon singlet solutions of Eqs. (^l|) for 
TV = 8 and < q < n. 
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singlets are more complicated than those of the triplets, 
the pattern is still recognizable to the extent that it may 
serve as a useful guide for applications to other system 
sizes. The reader will find it easy to identify in the tables 
of Part I one state for N = 4 and three states for N = 6 
which belong to this class of singlet excitations. In Table 
II we present the two-spinon singlet data for N = 8. 
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FIG. 5. Two-spinon (singlet) excitations with St = St = 
for iV = 16 (red circles). The configurations of (half- integer) 
Bethe quantum numbers (— | < I3 < . . . < Is < % ) are sym- 
bolized by the black circles and the Bethe quantum number 
1^ pertaining to the pair of complex solutions by black di- 
amonds. The N = 16 triplets from Fig. |lj are shown as blue 
crosses for comparison. 

Note that the pattern of {I3, . . . , /jv/2} is akin to the 
pattern of {Ii, . . . ,1^/2} f° r the triplets of a shorter 
chain. As N grows larger, the effect of the complex so- 
lutions z\ — Z2 on the energy (8a) relative to that of 
the real solutions Z3, . . . , z N / 2 diminishes and disappears 
in the limit N — > 00. The two-spinon singlets wilL.then 
form a continuum with the same boundaries @Jl3 Yet 
the effect of the complex solutions in the two-spinon sin- 
glets will remain strong for quantities inferred from the 
Bethe ansatz wave function (for example, selection rules 



and transition rates). 

How do we solve Eqs. (|3l|)? At q = tt the real roots 
23, . . . , z N / 2 of the two-spinon singlets are symmetrically 
distributed about z — 0. This observation paves the way 
for finding the complex roots exactly. Equation ( |31c| ) 
admits a solution in the limit |u| — » + ,|w| — > 1 + along 
a path with \u\ cx (\v\ — l) 1 '^. This limit matches the 
divergence on the left with that in the first term on the 
right. The Zj in the non-singular terms on the right have 
no effect on this solution. 

W hen we substitute the complex roots just found 
into (31b), the left-hand-side becomes ±Nir, while the 
sum on the right disappears bec ause of the symmetric 
Zj-configuration. Equatio n ( | 31b[) thus requires — 
(mod 2it). Finally, Eqs. ( |31a| ) with u = 0, \v\ = 1 can 
be solved iteratively similar to (|^), Rapidly converging 
solutions Z3, . . . ,Zn/2 are obtained for the configurations 
I3, . . . , In/2 indicated in Fig. [j| 

Significant computational challenges arise in the deter- 
mination of two-spinon singlet states at q ^ n, where the 
real roots are no longer symmetrically distributed and the 
complex roots have u > Q,u > 1. Now there is no way 
around solving Eqs. ( fil] ) simultaneously For this task 
a good root-finding algorithmic combined with carefully 
chosen starting values will be needed. 

The two-spinon triplets play an important role in the 
zero-temperature spin dynamics of quasi-lD antiferro- 
magnetic compounds. They are the elementary excita- 
tions of (|l|) which can be directly probed via inelastic 
neutron scattering. The two-spinon singlets, in contrast, 
cannot be excited directly from \A) by neutrons because 
of selection rules. The singlet excitations are important 
nevertheless, but in a different context. 

Several of the quasi-lD antifcrromagnetic compounds 
are susceptible to a spin-Peierls transition, which involves 
a lattice distortion accompanied by an exchange dimcr- 
ization. The operator which probes the dimer fluctu- 
ations in the ground state of ([!]) couples primarily to 
the two-spinon singlets and not at all to the two-spinon 
triplets. In a forthcoming column of this series, the focus 
will be on transition rates between the ground state of 
([!]) in zero and nonzero magnetic field and several classes 
of excited states that are important for one reason or 
another. 



V. SUGGESTED PROBLEMS FOR FURTHER 
STUDY 

1. Show that \Af±) is not an eigenvector of Ha- Cal- 
culate (Af±\S%S% +1 \J\f±) , a — x,y,x, to obtain the 
energy expectation value (M±\Ha\N'±) — —JN/4 
for the Neel state. Compare this result with the 
asymptotic ground state energy = (4 — 
\n2)JN inferred from (|l3|). 

2. Identify \A) for N = 4, 6 in Tables II and IV of 
Part I. Show that the state for N — 8 with Bethe 
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quantum numbers (||) and with real momenta 
fci = 1.522002, k 2 = 2.634831, k 3 = 3.64836, k A = 
4.761182 is a solution of (133) and (135). Then de- 
termine its energy E — Ep and wave number k. 

3. (a) Use the derivation of @ from (133) and (135) 
to establish the relation between the A^ and the 
Ii. Given that all A^ are integers, show that all 
Jj are integers for odd r and half-integers for even 
r. For the given (real) solution {z\ < ■ ■ ■ < z N / 2 } 
determine the one-to-one correspondence between 
the elements of (||) and (^). (b) Infer Q) from 
(130) with J replaced by —J and @ from (136). 

4. Calculate the {z^ of the state \ A) iteratively via (j|) 
from the {ij} given in ([?]) for various N. Compare 
the results for N — 8 to the {ki} given in Prob- 
lem |[ Check the energies for larger systems against 
the values listed in Table I. Extrapolate your finite- 
N data for (Ea — Ef)/ JN and compare your best 
prediction and its uncertainty with the exact result 

5. (a) Use (f|) and ( |l2[ ) to calculate the distribution 
of the magnon momenta in the ground state \A) 
via p (k) = J x oo dzcr (z)6(k - ir - 4>{z)). (b) Ex- 
press the energy per site of \A) as an integral of the 
magnon energy, — J(l — cosfc), weighted by po{k), 
and evaluate it to confirm the result (|l3|). 

6. In Problem || we interpreted \A) as a composite of 
N/2 magnons with momentum distribution (14). 
Increasing m z from zero means reducing the num- 
ber of magnons in the ground state to N(^ — m z ). 
From the solution {zi}, calculate the finite-TV dis- 
tribution pN(ki) = [N(ki — fci + i)] _1 of magnon mo- 
menta for m z = 0,0.125,0.25,0.375. Use data for 
N = 32 and for a much larger system, for example, 
N = 2048. Show that for h — 0, pN(ki) converges 
very rapidly toward (|l4]). Show that increasing h 
does not only deplete magnons, it also constricts 
the range of the magnon momenta. 

7. Gener ate data for E(S^) via (|). Calculate h(m z ) 
from ( |l9a| ) and Xzz(m z ) from ( 19b ). Then plot \ zz 



versus h for < h < 0.25 in relation to the exact 
value 1/(tt 2 J) at h = 0. Show that the slope of 
Xzz{h) increases without bounds as h — > 0. De- 
termine the Xzz at h = via extrapolation of the 
h 7^ data points, and compare your prediction 
with the exact result (EOT). 



8. The leading singularities in m z (h) and Xzz(h) at 
h = hg can be derived from the iV-dependence 
of the three highest energy levels in Fig. ||(a): 
E(S^)/J-N/4 = 0,-2,-4cos 2 (7r/2(iV- 1)), for 
Sf. = N/2, N/2 - 1, N/2 - 2, respectively. To calcu- 
late the third result, show first that z 2 = — z\ holds. 
Then derive the solution Zi = ±tan(7r/2(A r — 1)). 
Next use the relations h s = E(N/2) - E(N/2 - 1), 



h = E(N/2-l)-E{N/2-2),and m z = |-l/iVto 
express m z as a function of hs — h. The derivative 
of the result yields (||l]) . 

9. From Part I we know that the vector with Sf, = 
of a multiplet is related to the vector with — 1 
by an extra magnon with momentum k r = and 
an extra A,. = 0, while the {ki, . . . ,fc r -i} and 
{Ax, . . • , A r _ 1 } remain the same. According to (Q), 
this result implies z r = ±oo. Use (|^) to show that 
if the {zi, . . . , z r _i} are to remain the same for the 
two states, then the Jj-configuration must change 
as indicated in rows two and three of Fig. |^. 

10. Start from Eqs. (^) with Z\ = z* 2 = u + iv and 
real z$, . . . , zjv/2- Use 4>(x ± y) = arctan(2a;/(l — 
x 2 - V 2 )) + iatanh(2j//(l + x 2 + y 2 )) + 21-k with 
(undetermined) integer I, and employ the addition 
theorem 4>{zx) ± 4>{z 2 ) — 4>{[z\ ± ^]/[l ± z\z 2 ]) in 
sums and differences to separate real and imaginary 
parts of the resulting complex Bethe ansatz equa- 
tions. Take into account that I \ + 1 2 is always an 



integer. Adding N/2 to in ( |31b|) , which is per- 
mitted because of the undetermined 2Z7r's, makes 
1^ assume integer values near zero for the class of 
states of interest here. 
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